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Abstract 

We present a dynamical analysis of the (classical) spatially flat and negative 
curved Friedmann-Lameitre-Robertson- Walker (FLRW) universes evolving, (by 
assumption) close to the thermodynamic equilibrium, in presence of a particles 
creation process, described by means of a realiable phenomenological approach, 
based on the application to the comoving volume (i. e. spatial volume of unit 
comoving coordinates) of the theory for open thermodynamic systems. 

In particular we show how, since the particles creation phenomenon induces 
a negative pressure term, then the choice of a well-grounded ansatz for the time 
variation of the particles number, leads to a deep modification of the very early 
standard FLRW dynamics. 

More precisely for the considered FLRW models, we find (in addition to the 
limiting case of their standard behaviours) solutions corresponding to an early 
universe characterized respectively by an "eternal" inflationary- like birth and a 
spatial curvature dominated singularity. In both these cases the so-called horizon 
problem finds a natural solution. 



1 Introduction 



Though the "Standard Cosmological Model" [T], |[ finds its theoretical grounds 
on a purely classical and maximally symmetric solution of the Einstein equations, 
like the Friedmann-Lameitre-Robertson- Walker (FLRW) model, nevertheless it 
is commonly believed that quantum effects could have played an important dy- 
namical role during the evolution of the early universe. 

Among such quantum effects are surely of key relevance those related to pro- 
cesses of particles creation; indeed, due to the strong time variation singled out by 
the background gravitational field near the cosmological singularity, it is naturally 
expected that quantum fluctuations of the field itself, as well as the dynamics of 
any quantum field living on such a nonstationary background, lead to phenom- 
ena of matter creation, whose existence can in turn modifys the early universe 
behaviour. 

A detailed microscopic description of matter creation phenomena should be 
provided by a (fully self-consistent) quantum field theory on a curved space-time 
[[|, or, when the backreaction of the created matter becomes no longer negligible 
(in the sense of strong quantum fluctuations around any average effect), by a 
pure quantum gravity theory in presence of "matter" fields. 

A macroscopic phenomenological description of processes of matter creation, 
able to make account of the everage effect induced on the (classical) dynamics of 
an isotropic expanding universe, was proposed by I. Prigogine et al. (see 0-0) 
which is based on the application, to the universe comoving volume (i. e. spatial 
volume of unit comoving coordinates), of the standard theory of open thermody- 
namic systems having a conserved specific entropy (entropy per particle); such 
a description was then generalized in @ by allowing the specific entropy be no 
longer constant and putting its formulation in a completely covariant form. 

If the matter present in the universe is described by the energy-momentum 
tensor of a perfect fluid, then, in such an approach, the pressure of the fluid, apart 
from the usual thermostatic component, contains an additional term, directly 
related to the time variation of the particles number within the comoving volume; 
furthermore the total entropy within such a comoving volume, of the universe, 
directly proportional to the particles number, is clearly no longer constant (like 
in the Standard Cosmological Model), since it increases during the processes of 
matter creation (with interesting implications about the explanation of the actual 
universe entropy). 

Of course, the above mentioned macroscopic phenomenological approach is 
not completely self-consistent since it can not clearly determine the correct ex- 
pression for the rate of particles creation; therefore such an information is an open 
degree of freedom of the scheme and it should be provided by the formulation of 
an appropriate ansatz for the considered dynamical model. 

In H the influence of a particles creation term on a FLRW universe, having a 
flat three-dimensional space, was studied in presence of bulk viscosity, by showing 
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how the assumption of a reasonable ansatz (for the rate of particles creation) 
leads to the appearance of nonsingular cosmologies in correspondence to some 
particular choice of the free parameters of the model. 

Our analysis generalize the ansatz presented in || to make it even suitable 
for the description of a spatially negative curved FLRW model and the resulting 
early dynamics (analysed close to the thermodynamic equilibrium and in ab- 
sence of bulk viscosity, which, in agreement to what below discussed in section 2 
should be appropriately treated in this framework only as a small perturbative 
term) appears to be, in both the cases, as deeply modified by the processes of 
matter creation. Indeed in the case of a flat FLRW universe we get a nonsingular 
cosmology in correspondence to an open set of the allowed parameters. 

More precisely, by choosing an appropriate time variable (the logaritm of the 
comoving volume), we reduce the dynamical problem to an evolutive equation for 
the only energy density filling the universe, whose analytical solution provides a 
functional form of the time variable which is bounded or unbounded depending 
on the sign of the integration constant. In the physical case of a strictly positive 
integration constant, in correspondence to any choice of the allowed parameters, 
the evolution of the universe asymptotically approaches for sufficiently late times 
the FLRW one, while its early stages are characterized, respectively in the flat and 
negative curved cases, by an "eternal" inflationary-like behaviour and a spatial 
curvature dominated singularity. 

The most important aim of this paper is to show how, at least in the adopted 
phenomenological framework, the matter creation could have been deeply influ- 
enced the very early phases of the universe life, up to the point of providing 
an explanation for some fundamental puzzles of the modern cosmology, like the 
cosmological singularity and the horizon problem. 

It is important to observe how, though the precise nature of the phenomena 
of particles creation relays on the detailed microphysics working in the early 
universe, nevertheless the notion of heat Big-Bang (indeed we will get below 
the universe has a maximal initial temperature to be reasonably fitted by a value 
> 10 19 Ge^) implies that the most part of the created particles should be initially 
photons. 

In section 2 is given a brief review of the theory of open thermodynamic 
systems, whose application to the universe evolution provides the theoretical 
framework for the dynamical analysis developed in section 3; finally concluding 
remarks follows in section 4. 

2 Theory of the Open Thermodynamic Systems 

Here we provide a description of the theory of open thermodynamic systems 
and its application to the universe thermodynamics, which leads to the natural 
phenomenological description of matter creation processes adopted in the analysis 
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of section 3 0-0 (see also 0, §). 

Let us consider a thermodynamical fluid whose energy density, thermostatic 
pressure, volume, temperature, entropy and particles number are denoted respec- 
tively by p, p, V, T, S and N. All these thermodynamical parameters should 
satisfy the Gibbs relation: 

d(pV) = TdS - pdV + pdN ( 1 ) 

where \x denotes the chemical potential (i. e. the Gibbs function per particle) 
and is defined as: 

„ b + P) V TS (2 ) 
N N [Z) 

By the above definition of the chemical potential it is easy to rewrite the 
relation ([I]) in the form: 

TJ 

d(pV) = NTda - pdV + —dN (3) 

where a and H denote respectively the specific entropy (cr = S/N) and the 
total enthalpy (H = (p + p)V). 

Now, by the first law of thermodynamics, we get the fundamental relation 
(total energy conservation): 

d(pV) =5Q- PdV (4) 

where 5Q indicates the heat exchanged and P the total thermodynamical 
pressure. 

By the comparison of the two equations @ and (|) take place the following 
natural identifications: 

5Q = NTda P = p + Tl (5) 
where II denotes an additional pressure term due to the matter creation: 

„ H dN , s dlnN , , 

When applying this theory to the case of an expanding isotropic universe we 
should require that all the thermodynamic functions depend only on time and 
the volume V be the comoving volume corresponding (by convention) to an unit 
coordinate volume. 

Finally, by requiring that the universe be an adiabatic system, we get the 
following fundamental relation: 

SQ = NTda = (7) 
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Now, if we treate the universe as a closed system having a fixed number of 
particles per comoving volume (N = const.), then the above adiabatic condition 
(0 leads to the entropy conservation dS = (like in the Standard Cosmological 
Model); instead when the universe behaves as an open thermodynamic system 
and it is allowed the particles number per comoving volume changes during its 
evolution, the same adiabatic relation provides the following fundamental expres- 
sion for the entropy variation: 

rpn 

TdS = dN (8) 

To put the above thermodynamical scheme in a covariant form (appropriate to 
our purposes), it is enough to redefine the energy-momentum tensor of a perfect 
fluid as: 



T ik = (p + P)uiU k + Pg ik £,£ = 0,1,2,3 (9) 

where u\ and gi k denote respectively the four-velocity of the fluid and the 
space-time metric. 

To conclude, it is important to abserve that, since for dlnN/dlnV > (p + 
3p)/(3p + 3p) we have p + 3P < 0, then, in agreement to fli~l"l , the particles cre- 
ation term can provide a sufficiently negative pressure to remove the cosmological 
singularity. 

Let us now face the question concerning the possibility of including in the 
above scheme viscosity effects. 

Indeed it is commonly believed that phenomena of matter creation should 
be associated with the presence of a certain degree of bulk viscosity, as assumed 
in ||; however, as discussed in |J, when referred to rapidly space-time varying 
phenomena (like the ones concerning the very early cosmology) the standard 
theory of viscosity should be amended for accounting of the finite velocity by 
which any thermodynamical effect propagates. 

Such an aim can be reached, in a phenomenological approach, by introducing 
suitable relaxation times (typically of the order of the mean collision time) which 
make account for the characteristic time of the system evolution. When such 
characteristic time is sufficiently large, then the standard theory of viscous fluids 
can provide a satisfactory description of the considered thermodynamic system; 
but if such time becomes short enough, the system resembles more a viscoelastic 
solid and the causal effects should be taken into account. 

More precisely, when referred to the only presence of bulk viscosity, the theory 
of viscoelastic matter, derived in ||, can be (dynamically) summarized by the 
following energy-momentum tensor: 

T ik = (p + P + u)uiU k + (P + u)g ik £,£ = 0,1,2,3 (10) 
where u denotes the dilatational-stress coefficient and satisfys the equation: 
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LU + T U = ^ % ( ) = ( )X (11) 

with ; indicating the covariant derivative and r and £ respectively the relax- 
ation time and the dilatation coefficent; in the limit r —>■ oo the ratio £/r goes 
over the bulk modulus. 

Since such a theory is only appropriate to describe quasi-stationary systems, 
which evolve close to the thermodynamical equilibrium, then u should fulfills the 



inequality 10 



« P (12) 

Now the first law of the thermodynamics (which is provided by the conserva- 
tion law for the energy-momentum tensor ([H]) and in the below analysis plays a 
key dynamical role) reads: 



dp = -(p + P + uj)dlnV (13) 

Since u, r and £ are expected to be function of the only energy density p 
(below, by choosing a suitable ansatz for the particles creation term, we get that 
even P is a function of p), then a good (dimensional) assumption about the form 
the ratio £/r should take, is given by []: 

k 

p* = const. > k = const. > (14) 
Hence equation (|i~2f) rewrites: 

co « p 2 (15) 

By the dynamical analysis developed in the next section we will show that, in 
the case of physical interest, during the all universe evolution, the energy density 
satisfys the restriction p < p p = const. > 0. 

Now in the cases p <^ p* and p ~ p*, by (|i~5D, we see that in equation (|i~3|) 
the contribution due to to is to be regarded as a small perturbation. Instead in 
the case p ^> p*, the above claim is no longer always true. However the former 
cases appear preferrable from a physical point of view, since, in the limit of very 
small values of the viscosity coefficient £, we can yet speak of a real (weakly) 
viscous fluid (i. e. are yet allowed large values of the relaxation time r ). But the 
most important physical reason which leads us to believe the viscosity effects be 
negligible in the earliest phases of the universe evolution relays on the observation 
presented at the end of section 1 about the requirement of an early dominant 

x In [[To] it was assumed with satisfactory causal results the form £/tq ~ p. 
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radiation component; indeed termalized photons are not expected to resemble a 
viscous fluid, like instead could take place for a later hadronic component. 

According to the above general statement, in our approach we will neglect 
the viscosity effects which, providing only a small dynamical perturbation, can 
not drastically modify the results below obtained f\. 



3 Dynamical Model 

As well known, in a synchronous reference, the geometrical structure of the FLRW 
universe is summarized by a line element of the form : 

ds 2 = -c 2 dt 2 + R 2 (t)dl 2 (16) 

where c denotes the speed of the light, R(t) (t» < t < oo U > — oo) the 
cosmic scale factor (which has the dimension of a lenght) and dl 2 the line element 
of a constant curvature three-dimensional space, i.e.: 

dv^ 

dl 2 = , - + r 2 [dQ 2 + sin 2 9d(t) 2 ) (17) 

1 — kr 2 

Here by k we indicate the curvature signature {k = —1, 0, +1). 
The evolution of the FLRW universe is governed by the following two Fried- 
mann and "continuity" equations: 

2 (ldR\ 2 c 2 x c 2 k 

d( P R 3 ) = -Pd(R 3 ) (19) 

Above x = 8nG/c 4 (being G the Newton constant). 
Now we introduce the new time variable ri(t), defined by: 

dr] = 3H(t)dt R{r}) = i^e^"^ - oo < 77 < oo (20) 

where _R and r/ denote two assigned constants (R > 0). By other words 
we adopt as time variable, the logaritm of the comoving volume to describe 
the universe expansion (the above position dr] = 3H(t)dt clearly provides an 
appropriate time variable only in the case of an expanding universe, i. e. H > 0). 

In terms of (|20D the line element ([T6|) reads: 



<is2 = -^ + fl ° 2ef< "~' ,>)<if (21) 



2 However, on this level, we can not exclude at all that, under certain conditions, the viscosity 
effects could play a relevant dynamical role during the very early stages of the universe evolution. 
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Starting by the above line element and in agreement with the thermodynam- 
ical scheme presented in the previous section, the dynamical equations for the 
FLRW model, in presence of particles creation, take the form |7] : 



c 2 x c 2 k 



H 2 (V) = "-fp(v) - T^e-I^ (22) 

d T = -iP + P) (23) 
dn 

where p denotes the energy density associated with the matter filling the 
universe, and P the corresponding total (isotropic) pressure, constituted by the 
following two contributions: 

„ , ^ dlnN , , dlnN , . 

P = p + n n = - [p + P)li - m = - {p + P )— (24) 

Now, to the existence of such dissipative processes is associated an entropy 
production determined, in agreement with (H), by the relation: 

T dS = TSdN 
drj N drj 

Since, according to the second law of the thermodynamics, the entropy never 
decreases (dS > 0), from (p5|) we immediately get the conditions: 

dN > dlnN > (26) 
In terms of (PI ) equation (^) can be rewritten as follows: 

dp . ( dlnN \ . . 

By introducing the particles density n = -^e the above equation takes 
the more compact form: 

dp = p + p dn 
dr] n dr] 

The above thermodynamical scheme is completed by the following evolutive 
equation concerning the universe temperature 0: 

1_ d£ = fdp\ Idn 
T dn \dp J n n dn 

We conclude the formulation of the considered cosmological problem by ob- 
serving that the relation between the synchronous time and the adopted time 
variable rj is given by: 
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tM = t ° + Lsm ( » s(( "° ) (30) 

as well as, in terms of our variables, the physical horizon reads: 

d H ( v ) = ceh f V — ^ (31) 

The advantage in using the adopted time variable 77 consists of the decoupled 
character acquired by equation fl22|) (i. e. since now R(rj) is a known function, 
then equations ( p2|) amd (^7]) decouple from each other U), so reducing the 



considered cosmological problem to solve the only evolutive equation (|27j) . 

However to solve the equation ( |2~7|) it is necessary to assign the equation of 
state relating the pressure p to the energy density p, but overall to make an 
appropriate ansatz for the particles creation term dlnN/dt]. 

In what follows we shall assume an equation of state of the form [[13]: 



p=( 7 -l)p 1< 7 <2 (32) 

where 7 denotes the adiabatic index (in particular 7 takes the values 4/3 for 
the ultrarelativistic matter and 1 for a dust). 

In order to individualize an appropriate ansatz for dlnN/ dr], we preliminary 
observe that, by using (|32"D, then equation (^8|) immediately gives: 




(33) 

where po an d ^0 denote two positive constants (namely, po — p{Vo)- n o — 
It is easy to realize that equation (^) leads to the following implications: 

p = 0^0^^ = (34) 
drj 

dlnN 

p = const, fi =^> n = const, fi ^> — - — = 1 (35) 

dr\ 

Since in the case of a FLRW model with negative curvature (k = —1) the 
condition p = is yet compatible with the existence of a real space-time (char- 
acterized by H = cexp{ — (n — i] )/3}/Rq), we are naturally lead to the following 

ansatz : 

where a and (3 denote two positive constants and we set: 
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c 2 xct 

This expression for dlnN/dr] clearly satisfys 
p c = const., we should have: 

Pc = P 



and, in view of (|35|) , if p = 

(38) 



By using (^) and fl3"2|), equation ( p7[ ) can be rewritten in the form: 



f//9 



7P 



,P, 



1 



As it is easy to verify the solution of the above equation reads: 

p(v) = — r 

(1 + afrexp{'yP(r) - ?7o)}) /3 
or, equivalently, in terms of R as: 

p(R) 



3/37' 



(39) 



(40) 



(41) 



where a denotes an arbitrary integration constant. 

By assigning the initial condition p(r] ) = p , we immediately obtain: 



Po 



(1 + apP)* " " " ^ 

From equations Q33|) , ( f40l) and (^) we get the following expressions for the 
particles density as well as the particles number: 



n = n 



1 + af 



7/3 



1 + ap /3 exp{'y[3(ri - r] )} 



N = n 



1 + ap^expi'jP^rj - rj )} J 



(43) 



(44) 



Finally, by fl22|) , we get the following key expression for the expansion rate of 
the universe: 



c 2 x 



c 2 k 



\ 3 (1 + apf 3 exp{'yl3(ri - r?o)})^ R ° 



re 3 



(v-vo) 



(45) 



The expression ( |40|) shows how the energy density of the universe be a non- 
diverging function of 77 in correspondence to the choice a > (po < p), while for 
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a < (po > p) it diverges when 77 takes the value 77/ = 770 + [ln(— l/ap 13 )}/^^. 
However, for any value of the constant a, in the limit 77 — > —00, we get p — > p. 

As preliminary step, it is important to stress that the two cases a = (Bondi- 
Gold steady state universe) and a < 0, by describing solutions in which the energy 
density respectively remains constant or increases monotonically (diverging for 
77 = rjf), do not provide a suitable representation of our actual universe and 
therefore will be excluded from the following analysis. 

Then, it is easy to realize that in the limiting case when the processes of 
matter creation are absent (a = 0,i.e.p — > 00), we immediately reobtain, by 
(f45|), (|30|) and fl4l|), the standard FLRW dynamics. Furthermore our solution 
(a > 0) goes over the FLRW dynamics in the limit 77 — > 00, corresponding to the 
behaviour: 

P « ^ ( 46 ) 
Indeed the key difference between the Bondi-Gold steady state model and 
the our one, consists just of such asymptotic FLRW dynamics, which allows to 
overcome the two fundamental cosmological inconsistencies of a steady universe, 
concerning the existence of a microwaves background radiation and the unob- 
served level of particles creation, as required to keep constant the present energy 
density of the universe. 

More precisely, it is easy to realize that, by a suitable fmetuning of the free 
parameters p, (3 and a, it is possible to fit all the fundamental (predictive) features 
of the Standard Cosmological Model, including the nucleosynthesis process. 

Let us now analyse separately the two cases k = and k = — 1 corresponding 
to the choice a > 0. 

i) Flat FLRW model 

In the case of a flat FLRW model (k = 0) the function H(rj) behaves as the 
square root of p{rf) and, as it is easy to realize, it vanishes for rj — ► 00 while for 



rj — > —00 approaches the value cyxp/3. Thus, in the limit 77 — > —00, we get the 
following expressions: 



77 - 770 ~ cJ3 X p(t - t ) R(t) ~ R exp i c J*f(t-t Q ))~0 (47) 





p~p=^n~n = ri — iV~0 (4£ 



In addition it comes out easily from ( |30D and ([!!]) that the age of the universe 
and the physical horizon have a diverging character. 

Thus we see how the universe emerges from rj = —00 by an inflationary-like 
behaviour characterized by a (nearly) constant value of the energy and particles 
number densities, to be interpreted as the following asymptotic limits: 
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lb 11111 I "Tl/ I 

' T) — * — OO y Tj— * — OO \? 

being V the considered comoving volume (V = R 3 ) and U the total internal 
energy in it contained (U = pR 3 ). 

By other words, though U and A" vanish as t — > — oo, nevertheless their ratio 
to the vanishing comoving volume R 3 approaches constant values. 

ii) Negative curved FLRW model 

In the case of a negative curved FLRW model (k — — 1) we see from ( fi5| ) how 
the function #(77) vanishes for 77 — > 00 and diverges for 77 — > — 00. As it is easy 
to verify, for 7/^—00 take place the following relations: 



77-7/0-3/77 J; ^ ~ # + c(< - t ) ~ (50) 




N~Q (51) 



Above, the values p and n should be interpreted in terms of the limits ( (49| ) . 

By ( pl|) it is again easy to realize that the physical horizon diverges even in the 
case of a negative curved FLRW universe, while by fl3"U| ) the age of the universe 
is in this case finite ~ Ro/c + t — t . 

Since, asymptotically, the energy density approaches the constant value p, 
then we see how the very early dynamics results to be dominated by the spatial 
curvature (i. e. H <xl/R). 

In both the above cases (k = and k = —1), since, for 77 — > 00 (t — » 00) 
the energy density vanishes and the appropriate equation of state corresponds to 
7 = 1, then, in this limit, we get A^ ~ noRo 3 ^ 1 ^ / Po = NqoT 1 ^ / po. 

At the end we observe, how by fl25D , (|29|) and Cj32|) , the entropy and the 
temperature of the universe can be expressed respectively in terms of the particles 
number and number density as follows: 

S = S ^- T = tJ-Y 1 (52) 



N \n 

By the above expressions for N(rj) and 77(77), it is easy to construct the explicit 
functional form for S(rj) and T(r)). In particular it is worth noting how in the 
two cases k = and k — — 1 the entropy reaches for t —>■ 00 the (maximum) value 
S ~ Sqcl^ 1 ^ I po, while the (maximum) initial value taken by the temperature is 
given by: 




Tmax —Tq\ — (53) 
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In particular for an initial radiation dominated universe (7 = 4/3) we get the 
maximum value T' max = T (p/ po) 1 ^ 4 - 

Of course the values of the above maximal temperature can be suitably fixed 
to avoid any contraddiction with the observed universe. 

4 Concluding Remarks 

Finally we enphasize how the analysis of a positive curved FLRW universe (k = 1) 
can not be developed in terms of our approach since the time variable r] would 
be hill-defined in the "turning point" of the universe (H = 0) (as well as not 
applicable to the collapsing picture H < 0), naturally expected in the evolution 
of such a model. Furthermore, as it is easy to realize, for negative values of H 
the adopted ansatz would lead to a negative synchronous time rate of particles 
and entropy production, contraddicting the second law of thermodynamics; for 
these reasons, the analysis of this case requires further investigation. 

As the results of our analysis show, the presence of a particles creation term 
deeply influences the universe dynamics, by allowing the energy density and the 
temperature have never diverging character. More precisely we get, as above 
outlined, a drastic departure from the "Standard Cosmology" only in the (very) 
early stages of the evolution, while for sufficiently later times (depending on the 
adopted choice of the free parameters), the obtained dynamics approaches the 
usual FLRW one. 

It is worth enphasizing how, as a direct consequences of this departure, the 
obtained solutions for the two cases k = and k = — 1 both provide, by the 
diverging character of the physical horizon, a solution to the so-called horizon 
problem. 

We conclude by observing how any physical motivation for the adopted ansatz 
is rather difficult to be provided since it should be referred to an average pro- 
cedure over all quantum effects being active in the early universe. In particular 
the impossibility to make precise account for the quantum gravity effect lim- 
its strongly the physical derivation of the particles creation rate (indeed, in the 
framework of our phenomenological approach, the matter is created by the gravi- 
tational field itself by its time variation and the main contribution to this average 
effect can be expected to come from the quantum fluctuations of the cosmological 
background) . 

However, in spite of this great difficulties, the relevant influence the matter 
creation phenomena could have exerted on the very early universe dynamics, we 
believe, makes of interest any effort in this direction. 

Daniele Oriti and one of the two referees are thanked for their valuable advices 
on this matter. 
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